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We consider the conformal Einstein equations for massless collision- 
less gas cosmologies which admit an isotropic singularity. It is shown 
that the Cauchy problem for these equations is well-posed with data 
' consisting of the limiting particle distribution function at the singular- 
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1 Introduction 



In the first two papers of this series (Anguige and Tod 1998a, b, henceforth 
ATI, ATII) consideration was given to the Cauchy problem for the confor- 
mal Einstein equations near an isotropic singularity for two different matter 
models: the perfect fluid (ATI) and the collisionless gas (ATII). This prob- 
lem was solved in full generality for the perfect fluid but for the collisionless 
gas we had to assume spatial homogeneity to make progress. In this paper 
we again take the collisionless gas as matter model, the aim being to extend 
the existence and uniqueness results of (ATII) to the inhomogeneous case. 

Recall that a spacetime (M,g a b) has an isotropic singularity if there 
exists a manifold M D M, a smooth Lorentz metric g ab on M, and a func- 
tion Z defined on M such that 

g ab = Z 2 g ab for Z > (1) 
Z -> on S (2) 
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for some spacelike hypersurface S in M. Given some matter model one seeks 
to generate cosmologies with isotropic singularities by solving the Cauchy 
problem for the conformal Einstein equations, with data given for the un- 
physical metric g a b and the matter at the singularity S. 

In (ATI) it was shown that with a polytropic perfect fluid as matter 
source the free data for the conformal field equations consist of just the 3- 
metric of the singularity: the second fundamental form of the singularity 
has to vanish and there is no data for the matter. Every choice of 3-metric 
then determines a unique 7— law cosmology with an isotropic singularity. 

With a collisionless gas as source the field equations of general relativity 
are known as the Einstein- Vlasov (EV) equations, for the metric g a b and the 
particle distribution function /. In (ATII) it was shown that the massless 
EV system transforms nicely under a conformal rescaling as in (l)-(2) and 
that one may identify the free data for the field equations at the singularity 
surface S. The situation is rather different from the fluid case. One is free 
to prescribe the limiting particle distribution function f° at the singularity, 
subject to a single integral constraint. This function then determines the 
first and second fundamental forms of the singularity surface. The second 
fundamental form need not be zero here. 

Under the assumption of spatial homogeneity for the metric and the mat- 
ter it was proved in (ATII) that each f° determines a unique EV cosmology 
with an isotropic singularity. In the present paper we show that the same 
result holds in the absence of any symmetries. 

With a convenient choice of conformal gauge and in harmonic coordi- 
nates, the inhomogeneous field equations may be written as a symmetric 
hyperbolic system of integro-differential equations with a ^ forcing term, 
for certain carefully chosen matter and metric variables. These equations 
are essentially of the form studied in (Claudel and Newman 1998), but the 
matter integrals present in the Einstein equations prevent a direct applica- 
tion of the existence and uniqueness result obtained there. It is however 
possible to combine certain techniques used in this work with the standard 
method of energy estimates to obtain the desired result. We note that the 
semi-group theory used by Claudel and Newman is entirely avoided in our 
approach. 

The paper is organised as follows: 

In section 2 a brief review of the Einstein- Vlasov sytem and its behaviour 
under conformal rescaling is given. 

In section 3 we collect a few results from (ATII) on the initial data for 
the field equations and on conformal gauge fixing. 

Section 4 contains a proof of the main result, which is summarised in 



2 



Theorem 4.1 at the end of the paper. 



2 The massless Einstein-Vlasov system and con- 
formal rescaling 

A collisionless gas in GR is described by a positive function f(x a ,pb) on 
the spacetime cotangent bundle, which represents the number of particles 
near the point x a with 4-momentum near p b (Ehlers 1971). The condition 
that the matter be collision-free is equivalent to requiring that / be constant 
along the geodesic flow, which is the vector field defined, in local coordinates, 
by 

ab d l dg ab d 

The statement Cf = is known as the Vlasov equation. 

For massless particles the function / is supported on the seven-dimensional 
submanifold of the cotangent bundle given by the equation g ab p a Pb = 0. The 
stress-energy-momentum tensor due to massless particles is given by 



<ab = / fPaPb^- Q d 3 p (4) 



with the positive quantity p° being determined by the relation g ab p a Pb = 0. 
Note that the massless condition implies T a a = 0. 

If the Vlasov equation holds then this stress tensor is divergence-free 
independently of the Einstein equations being satisfied. 

The coupled Einstein-Vlasov equations, for the metric g a b and the par- 
ticle distribution / are 

G ab = 8tt / fpaPb^-^ d 3 p (5) 



CgU) = (6) 

Suppose now that we have a physical spacetime (M,g a t,) and a massless 
particle distribution function / so that g a b and / satisfy tilded versions of 
(5)- (6). Suppose also that we have an unphysical spacetime (M, g a b) defined 
by 

g ab = Z 2 g ab (7) 
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for some conformal factor Z on M D M. 

As discussed in (AT) it is a consequence of the conformal invariance of 
null geodesies that the massless Vlasov equation is also conformally invari- 
ant. This is to say that if we simply write / = / then the Vlasov equation 
for / in M implies the Vlasov equation for / in M: 

Cgf = (8) 

If we now use the unphysical metric g ab to define an unphysical stress tensor 
via 



then 




and one may now use the unphysical Vlasov equation (8) to show that 
V a T ab = in M, just as V a f ab = in M . 

The conformal EV equations for g ab and / are now 

R ab = 2V a V 6 log Z - 2V a log ZV 6 log Z 



+ g ab {D log Z + 2V C log ZV C log Z) + -^T ab (11) 

and 

£g(f) = (12) 

From now on we assume that a rescaling of the spacetime metric as in (l)-(2) 
exists, and our aim will be to solve equations (11)-(12) for / and g ab in the 
neighbourhood of the singularity surface X. 

In (ATII) it was noted that the massless condition on the matter implies 
that the conformal factor Z must remain smooth at S. In what follows we 
will require Z to be a good coordinate in M, and we thus make the following 
extra assumption: 

Assumption 2.1 The conformal factor Z is such that V a Z / at E. 
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3 Initial data and conformal gauge fixing 



The singular nature of the conformal field equations (11), (12) imposes con- 
straints on the data at an isotropic singularity S ,which are rather different 
from the usual constraints on a regular hypersurface. To say what these are 
we first make the following definitions: 

Let h®j be the 3-metric of the singularity surface S and let Kfj be the sec- 
ond fundamental form of £ in M. f° will stand for the orthogonal projection 
of / onto £ : fix^p^) = f{0,x\p b ) and V° = [g ah V a ZV h Zfl 2 , as eval- 
uated at S. 

The constraints then read as follows: 

f°(x i , Pj ) Pk d 3 p = (13) 



^-^J wrZ.)* *' (14) 



where K° = (h^Kfj , 

4vr f fpiVjPkPi , 3 



Xijkl = nvP J ((/.rw 2 p ( j 



4?r /" fpiPjPk , 3 / ir , 
(yO)Wi ((^PmPr.) P ^ j 

and -Dj is the covariant derivative operator associated with h^. 

Now (13) is just an integral constraint on f° and from (ATII) we have the 

following theorem: 

Theorem 3.1 Let f°(x l ,pj) be a smooth, positive function on U x M 3 , 
where U is an open subset of M 3 . Suppose that for each x G U 

1. f° is compactly supported in p. 

2. /° is supported outside some open ball containing p = 0. 

3. /° is not identically zero in p. 
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Then, given a smooth strictly positive function V on U, there exists a 
unique, smooth, positive definite 3-metric h - on U satisfying (14) and a 
unique smooth trace-free symmetric tensor Kfj satisfying (15) 

Also from (ATII) we have that the freedom in the choice of unphysical 
metric g a b may be used to obtain the following: 

Theorem 3.2 If (M,g a b) is a solution of the massless Einstein- Vlasov 
equations, with an isotropic singularity, and Assumption 2.1 holds, then the 
conformal factor Z may be chosen so that V° = 1, K° = 0, and Z is a 
harmonic function in M. 

We will henceforth refer to the gauge choice of Theorem 3.2 as the 
harmonic gauge. 

In summary the initial data set for the conformal EV equations at Z = 0, 
in the harmonic gauge, is as follows: 

• f° is free data, subject only to the integral constraint (13). 

• The initial 3-metric is determined by f° via (14) 
with V° = 1. 

• The initial second fundamental form Kfj is determined by f° via (15) 
with V° = 1, K° = 

4 The Cauchy problem for the conformal EV equa- 
tions near S 

In this section we show that the harmonic gauge conformal EV equations 
may be solved by a combination of the standard method of energy estimates 
with certain of the techniques used in the proof of the Newman-Claudel 
theorem. This theorem cannot be applied directly because of the presence 
of matter integrals in the field equations. 

4.1 The reduced Einstein equations as a symmetric hyper- 
bolic system 

The first step in solving the field equations is to impose the harmonic coor- 
dinate condition which leads to a system of reduced equations. 
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Given a coordinate system we make the following definition 

Rpv = Rfiv + 9a(n9 v )H a (18) 

where H a = \Dx a . 

The reduced conformal EV equations are then defined to be 

Rjlu - |v M v,z + A Vm zv,z - ^„ = ir,, (19) 

c g f = (20) 

If we put x° = Z and /i a/ 3 7 = g a /3^, then equations (19)-(20) can be decom- 
posed as 

d z g a p = h a (30 (21) 
- 9 ij d z h aPj = -g ij djh a p (22) 

dzg a(3 = -g^g^ h ^ (23) 

g 00 d z h afS0 = -g lJ dih a/3j - 2g 0l d i h a p ld + F a p(g,g- 1 ,h) 

+ + 2<7° 7 (/l QT/3 + hfra - K Pl )} (24) 



<'°% + ^-h a ^^wr° (25) 



where 

V = |{4V M ZV,Z - <?°V - T M 4 (26) 
and F a p is polynomial in its arguments. 

It follows that the reduced Einstein equations (21)-(24) are, for fixed /, in 
symmetric hyperbolic form. 
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4.2 Initial data for the reduced equations 

We would like to solve the conformal EV equations in the harmonic gauge. 
So, as initial data for the reduced equations first choose a positive smooth 
function /° on T*£ satisfying the constraint (13). We also require f°(x i , •) 
to have fixed compact support and to be supported away from a fixed neigh- 
bourhood of the origin for x % in a compact set. Then put (g ) 00 = (V ) 2 = 
1 and choose = 0, so that the x % are 'initially comoving'. We may 
now use Theorem 3.1 to determine gfj and h^ = 2Kfj from (14) and (15) 
with K° set to zero. 

From (19) one must have that \3Z = initially. But by definition 

DZ = d a g a0 + \g a0 g p(T d a g p<J (27) 

so that hooo = initially. 
Equation (24) implies 

Z a/3 = 2 { /l a/30 ~ h a0f3 ~ h % a } ( 28 ) 

so that Z% = 2h° ij0 , Z$ = Zl = 0. 

In order to recover the full harmonic gauge equations from the reduced 
equations we will have to show that H a = \Dx a = in M. We therefore 
pick /i? 00 = and choose spatial coordinates x % satisfying the following 
condition on S: 

g mn r mn = (29) 

where Y l - k are the Christoffel symbols of g®j. These choices ensure that Dx l = 
on E. 

4.3 Recovering the harmonic gauge equations 

We now use the contracted Bianchi identities to show that if {f,g a p) is a 
solution of the reduced equations (19)-(20) with data as prescribed in 4.2, 
then DZ = Dx l = and (f,g a p) solves the conformal EV equations (11)- 
(12). 

First write 

-V = TP2^ V + |v M V,Z + ^„ - A V M ZV,Z (30) 
We know from (20) that V^T^ = and this leads to 

VS^ = ^(V a Z)(R au - S av ) + ^d u H° - ^(d u Z)H° (31) 
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and now (19) gives 

VUj, = ^{V a Z)(-gp {a d u) HP) + ^d„H° - ^(d„Z)H° (32) 



From (19) one has R H = §#° , so that 



R=^H°-g^g l{a d /3) H^ (33) 



The contracted Bianchi identities, V^Guj, = 0, now imply 



= (V a Z)(-g p{a d u) HP) + -d u H° - — 2 {d v Z)H« 



- V a {9 P{ ad v )HP) - \v v (^H° - g^g l{a d p) HA + F u (g, g' 1 , h, H, dH) 

i \z j (34) 

where F v is polynomial in its arguments and homogeneous in H, dH. 
Now multiply (34) by g vX to get 

= ~\sT^H x - hv a Z)g^g p{a d v) HP 



+ ^g vX d u H° - lg" x (d v Z)P° + F x (35) 

where P° = ±H°. 

From (35) one gets the following at Z = 0: 

- 2(V a Z)g» x gp {a d u) HP + g vX d u H° - 3g» x (d v Z)P° = (36) 

Now P°(0) = d z H°(0), so (36) implies that P° = d z H° = at Z = 0. 
Equation (36) also implies that d z H l = at Z = 0. 
From the definition of P° it follows that 

If we put /ig = 3 a i^ then (35) and (37) can be written in first order form 



as 



■ dhf ..dhS 

s "nk'-o'ni (38) 
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(39) 



(g 00 )^ = 2g°%h x -g ij d i h x 

-^{g° a h X a + Sg 0X P°} + F x (41) 

The system (38)-(41) can be written in an obvious way as: 

a\u)—=a\u)— i +b(u)u + -c{u)u (42) 

with a +ve definite and a a symmetric. 

The eigenvalues A of ((a°) _1 c)(0) satisfy either 

A = (43) 

or 

AP° = _p° + (44) 

\h% = -2g 0l3 (0)h%-6g 0a (0)P° (45) 
Putting a = i in (45) gives 

A/4 = -2/4 

while putting a = in (45) gives 

A/4 = -2/4 - 6P° (46) 

So if A / 0, -2 then 

(2 + A)(l + A)+6 = (47) 

Thus (a°)^ 1 c(0) has no positive integer eigenvalues. Now, as remarked in 
section 4.2, the choice of initial data for the reduced equations forces the 
quantities \Dx a to vanish at Z = 0. Therefore (42) has the trivial solu- 
tion u = and by Theorem A.l (see Appendix) this is the only smooth 
solution with the given data. It follows that any smooth solution of the re- 
duced equations (19)-(20), with data as in 4.2, is a solution of the harmonic 
gauge EV equations . 
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4.4 Solving the reduced field equations 

We aim to solve the reduced EV equations on a manifold of the form M x 
[0, T], with M a paracompact 3-manifold. We will use the method of energy 
estimates to obtain local (in space and time) solutions of the field equations, 
which can then be patched together using the finite domain of dependence 
to obtain a solution in a neighbourhood of the singularity S = M x {0}. 

4.4.1 Localisation 

First choose a locally finite cover of £ by relatively compact harmonic co- 
ordinate charts O a . Now construct a second cover {0' a } in such a way 
that 0' a C O a . 

Restrict attention to a single chart O a and consider the reduced equa- 
tions (19)-(20) on O a x [0, T]. In order to obtain energy estimates it is 
desirable to have data in C^°(0 Q ). We thus cut off /° to zero smoothly 
along with g^ + Sij. We now work with g a p = g a p — r] a p as field variable 
{Vafi = diag(l, —1, —1, —1)). Obtain Kfj from the new /°, g a p via equation 
(15). This cutting-off procedure is done in such a way that the original data 
is maintained on D 0' a and is zero outside O'" D 0'^. 

Now note that the singular field equation (26) will no longer be consistent 
at Z = with the localised data. To remedy this we introduce a 
function F(x i ) € C£°(O a ) such that < F(x i ) < 1 and 

F(x) = l° xe °*\°° - (48) 

We use the function F to obtain the following localised modification of (26)- 
(27): 

R% = Ftf) \\^^vZ + -^(T^ + g 00 g^ - 4V M ZV,Z)1 

{Z Z j (4g) 

L g f = (50) 

We will solve (49)-(50) on O a x [0, T a \ and then use the finite domain of 
dependence to argue that this gives rise to a solution of conformal Einstein- 
Vlasov on 0' a x [0, T a ], determined by the right data. 
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4.4.2 Approximate solutions of the field equations 

An important step in solving the reduced field equations is, as in the proof 
of the Newman-Claudel theorem, to write them in such a way that they are 
manifestly formally solvable. One then obtains approximate finite Taylor 
series solutions which solve the equations up to 0(Z q ), q arbitrary. We note 
that the hyperbolicity of the field equations is not used at this stage. 

First define 

V = 4 {4V„ZV„Z - 5 °V " - (4V^ZV U Z - 5 °V " ^)(0)} 
Z (51) 

Then (49) is just 

R% = *V) {|v M v,z + ± v} (52) 

(F was chosen to be zero wherever (4V ^ZV U Z — g m g^ u — r^)(0) / 0) 
Define also 

Z a p = Z a p - Z a p(0) h a/ 3j = h af 3 y - /i Q/ 3 7 (0) (53) 
In terms of these variables, equations (49)-(50) become 

- g ij d z h al3j = -g ij dj(h a pQ + h aP0 (0)) (54) 

d z g a p = h a/ 3o + h a/ 3o(0) (55) 

d z S afi = -9 Pv 9°" i {h ltv0 + Vo) (56) 

g°°d z h a(S0 = -g ij di(h a/3j + h af3j (0)) - 2g 0i d i (h afS0 + h a/30 (0)) 
+2F a p + 2F(x i ){F 1 y(h a ^ + h Pia - h afh ) 

+ ^p-{Z a p + (HtfWihcryf, + h Ma - h aPl )} (57) 

and 

L g f = (58) 
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where 

(F 1 y(Z) = 1\{H{P)'{SZ) dS , {H{p = g<* (59) 
Jo 

The prime denotes differentiation with respect to time, and F\ can be com- 
puted from (53) and (56). 

Meanwhile we calculate the following evolution for 

8 Z Z^ = -F 2 (Vo + Vo(0)) + (F 3 ) a ^(hapo + h af so(0)) 
+(^4)^(^0 + h apo (0)) + (F 5 )^ + (F 6 ) a ^(h a/3i + h a pi) 

-{F7) a ^ v {h a i3Q + h af 3 ) 

+|{ - v - i? 2 (0)Vo + (H 3 ) a ^(0)Kp 

+ (^4)^(0)^/30 + (^6)^(0)^ - (^ 7 )^(0)^o } (60) 

where Fj is to iZj as Fi is to Hi and 

#2 = 2g 00 (61) 

(ir 3 )°£, = 2 M °V' 3 (62) 
mt, - (det^)- 1 / 2 ^ / lM± <i 3 p (63) 

(H 5 ^ = 2 { det 9 )-^J^§t id ^ (64) 

{H^^ { ^ g) -^J P J^^ P (65) 

(# 7 )f„ = (det^- 1 / 2 / ^( p V((P0)-W+P^+^°)) <* 3 p 

J {P j (66) 
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The Fi are regular terms, and can be calculated from (54), (55), (58). 
The whole system of equations (54)- (58), (60) can be written in the form 

A°(u)d z u = A l {u)d iU + B(u)u + -=C(x\ u)u (67) 

(i=l,2 ... 6) for some matrices A a ,B,C with A positive definite. The 
vector u stands for {f,g a p,g aP ,h a ^,Z a p). 

Lemma 4.1 The matrix C(x l ,u) in (74) satisfies C(x l ,u)u° = Vu and 
(A°)~ 1 C(x t , u ) has no positive integer eigenvalues. 

Proof Since Z a p(0) = h a ^(0) = the first part follows. Now note that 
the eigenvalues of ((^4°) _1 C)(0) satisfy either A = or 

A/i Q/ 3o = F{x l ) {Z af5 + {h a0f3 + hf3 0a - h af 3 }} (68) 
XZ afS = -Z a p + 2 j - K po + 9 a fs(0)g^(0)g 0u (0) Vo 

+ 1 -(detg(0))-^g^(0)h, u0 j d 3 p 

1 J [P > J (69) 

where p = p° = (~g lj (0)PiPj) 1/2 and h af3i = 0. 

Now if F(x l ) = then either A = or /i Q/ go = 0- If A / then from (68)- 
(69) we must have A = — 1. It is therefore enough to consider those x at 
which F(x) / 0. At such x the relation between /° and gij{0) is given by 
(14). 

Equations (68)-(69) can be expanded as follows: 

A/iqoo = F{Z 00 + 2h 000 ) (70) 



AZ 00 = -Z 00 - 6/iqoo + 6g lj (fyhijo 



(71) 



A/ijoo — F(Zu 



XZn 



= -Z l0 - 2h m - 2(detg(0))- 1 / 2 h jk0 J j^^Pi d 3 P 



(72) 



(73) 
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Xh ij0 = F{Zij - 2h ij0 ) (74) 
XZij = -Zij - 2h ij0 + 2^(0)/i 00 o - gij(0){h mn0 g mn {0)) 

-^nno(xT) ( 75 ) 

where 

X "™ = (detff(0))- 1 /2 J _£_ p rn p n p . p . ^ p (?6) 

Taking the trace of (74)-(75) leads to 

(0 ron (O)/w>)((A + 1)(A + 2) + 6F) = 6Fh 000 (77) 
while (70)-(71) give 

((A + 1)(A - 2F) + 6F)/i 00 o = 6Fg mn (0)h mn0 (78) 

and hence 

(g mn (0)h mn0 )({(\ + 1)(A - 2F) + 6F)((A + 1)(A + 2) + 6F) - 36F 2 ) = 

(79) 

Now if A were a positive integer then the polynomial in F and A on the 
left hand side of (79) would be greater than or equal to 12(1 + 2F(1 — F)). 
Since < F < 1 it follows that if g mn (0)h mn o / then A is not a positive 
integer. 

Suppose now that g mn (0)h mn o = 0. Equations (74)- (75) then become 

\h ij0 = F{Zij - 2h ij0 ) (80) 

XZ^ = —Zij — 2hij + 2<7jj (O)foooo - ^mnO^X^j 1 ) (81) 
and (70)-(71) become 

Xhooo = F(Zqq + 2/iooo) (82) 

XZqq = —Zoo — 6/iooo (83) 

which imply 

((A + l)(A-2) + 6F)/i 000 = (84) 
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Suppose (A + 1))(A - 2) + 6F + 0, so that h 000 = 0. Then (74)-(75) imply 

((A + 1)(A + 2) + 2F)h ijQ = -Fh mn0X ™ (85) 

But clearly the eigenvalues of x are positive and hence A cannot be a positive 
integer. 

Suppose finally that h mn o = 0. Then (72)-(73) are 

Xh m = F{Z i0 ) (86) 

XZio = —Zio — 2/tjoo (87) 

which imply 

(A 2 + A + 2F)h m = (88) 

and we conclude from all this that ((A°)~ 1 C)(0) has no positive integer 
eigenvalues. 

Lemma 4.1 now allows us to construct approximate solutions w of the 
field equations (67). Specifically we make the following ansatz 

g 7P 

w = Y^ —™ {P) with w(0) = w {0) = u° (89) 
p=0 p ' 

and now note that 

Z(A°(w)d z w - A i (w)d i w - B{w)w) - C(w)w = 

q ZP 

—(pA°{w) - C{w))w {p) - Z{A\w)d iW + B(w)w) (90) 

Let the be determined by demanding that the Z-derivatives up to or- 
der q of the right hand side of (90) vanish at Z=0. To see that this prescrip- 
tion is well-defined note for example that the first derivative is 

1 vp—l 1 yp 

y -=- — ttMV) - c( w ))w^ + y ±-{pA Q [w) - c(w)Yw^ 

z — ' (v — 1 ) pi 

p=l yF ' p=0 F 

-Z(A i (w)d i w + B(w)w)' - (A i (w)d i w + B(w)w) 
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Setting Z = and equating the resulting expression to zero gives 

(A°(u°) - C(u°))wW - (C{w))'{0)u° - (A^u^diU + B(u°)u°) = 

Lemma 4.1 now implies that is uniquely determined by the initial 
data u°. Equating the pth derivative of (90) to zero at Z = gives an 
equation of the form 

(pA°(u°) - C(u°))w ip) + 4>{w {p ' l \ . . . ) = 

Thus the are inductively well-defined and are Cfi° by the choice of initial 
data. 

By Taylor's theorem we must now have that 

Z(A°(w)d z w - A i (w)d i w - B(w)w) - C(w)w = ( (91) 
for some ( € admitting an estimate ||£|| s < K s Z q+1 , say. It follows that 

A°(w)d z w - A i (w)d i w - B(w)w - ^C(w)w = £ (92) 

Z 

where £ 6 for each Z and admits an estimate ||£|| s < n s Z q . 
We now wish to check that the approximate solution w = (f q , g g , g~ l , dg q , Z q ) 
of equations (67) gives rise to an approximate solution of the original lo- 
calised reduced equations (56)-(57). First a calculation shows that the pro- 
cess of writing 

Z-^H^Z) - Hi(0)) = f 1 H' l (SZ)dS 

Jo 

which led to (67) can be inverted to obtain an approximate solution of 
equations (50) and (52) with error 0(Z q ). Now one needs to check that the 
part Z a p of the approximate solution is close to the original definition (51). 
To do this define Z^ by the right hand side of (51) expressed in terms of the 
appropriate parts of our existing approximate solution. Then one calculates 
that 

d z {z^ - v(o)) = -|( v - + §{* + °( zq )} ( 93 ) 

where * is most of the right hand side of (60): 

d z (Z, v - V(0)) = " V(0)) + §{*} + 0(Z q ) (94) 

Subtracting (93) from (94) gives 

d z (Z(Z iiV -Z ilu )) = 0{Z q ) (95) 
Hence Z^ - Z^ = 0(Z q ). 
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4.4.3 The reduced EV equations as a coupled symmetric hyper- 
bolic system 

While the form (61)-(65), (67) of the reduced equations is useful for demon- 
strating formal solvability, it is not appropriate for the application of the 
method of energy estimates. The presence of the terms in the equation 
for Z^y spoils the symmetry of the system as a whole. On the other hand, 
while the field equations in the form (49)-(50) are symmetric, the terms 
are too strong to be handled by the contraction mapping technique of 
Claudel and Newman. We need to choose variables in such a way that 
the field equations are symmetric as a whole and have singular terms which 
are no worse than This can be achieved by defining the following new 
variable: 

P a p = </>(W a ZVpZ - g m g a p) - ^-p fp aPf3 (96) 

where <p = (j)(pi) G Cq° and J R3 4> d 3 p = 1. We also suppose that supp</> C 
supp/°(x, •). Now define 

Qaf3 = ^(PafS ~ P a( s(0)) (97) 

with Qa/3(0) = 2d z P a p{Q) determined by the localised data. In particu- 
lar Q^(0) e C °°(O a x R 3 ) and suppQ^(0,^, •) C suppf°(x\ •). 
The reduced field equations now decompose as 

d z g a (3 = ha/30 (98) 
- 9 ij d z h aPj = -g ij djh a p (99) 

d z g a(S = -ffVVo (100) 
g 00 d z h a(30 = -g lJ djh a pi - 2g )l d i h am + F a(3 (g, h) 

F(x i 



+ 



2g 0l (h aJ f3 + hp la - h aM ) + I Q a /3 d 3 p\ (101) 

Jm. 3 ) 
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Using the Vlasov equation one calculates the following evolution for Q a p 

dzQa "--^^ + 2 {dt9 ] ^-dp7 

- 2^d t P a ,(0) - »(^)|- (^defcV) (0)1 

P° P° dpi V P° /J ^ 103 ^ 

Clearly the system (98)-(103) is symmetric hyperbolic. 

We now note that the terms containing p° in equations (102)- (103) blow 
up at the vertex of the light cone and thus cause technical problems there. 
For this reason we multiply the right hand sides of these equations by a 
function ip(pi) 6 Cq°(IR 3 ) which is equal to unity on a set containing the 
support of f° and equal to zero in a neighbourhood of the origin. Call the 
resulting equations (102)', (103)'. 

We now wish to show that the approximate solution (f q ,g q ,h q ) of the 
field equations obtained earlier gives rise to an approximate solution 
(f q ,g q ,h q ,Q q ) of (98)-(101), (102)' - (103)'. To this end define Q a/3 = 
2Z~ l (P a p — P a p(Q)) with P a/ 3 constructed from f q ,g q via equation (96). 
If (f q ,g q ,h q ) satisfies (56)-(57) up to 0{Z q ) then one has (dropping the 
subscript q for the moment) 

g 00 d z h a/3 o = -g %3 dihapj - 2g 0t dih a/30 + F a/3 
{ [ Qa? d 3 p + 2g^< (h al p + hp la - hafr) ) + 0{Z«) 

J (104 ) 

A calculation gives that this Q a/ 3 satifies 

d z Q a p = rhs + 0(Z 9_1 ) (105) 

where 'r/is' is the right hand side of equation (141)'. It follows that we have 
an approximate solution (f q -i,g q -i,h q -i,Q q -i) of the field equations with 
error 0(Z« _1 ). 
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4.4.4 The field equations in matrix form 

Let the vector u = u(Z,x % ) stand for (g a p, g a/3 , /i a/ 3 7 ). Then we may write 
(98)-(101) in the following form 

Q 



i=l K 



//)+ / | ) <f'l> 



(106) 



where A a ,G,C are polynomial in the components of u, A a are symmetric 
matrices and A is positive definite. 

Let v = v(Z, x l ,pj) stand for (/, Q a p). Then we may write (102)' — (103)' 

as 



6 1 
= Y, Bl (Pj^)d l v + -D(u,v, Pj ) (107) 



dv 



dZ 



where the B l are symmetric matrices and B l , D depend smoothly on Pj,v and 
on u away from detg = 0. 

Now let u q and v q stand for the appropriate parts of the approximate 
solution we obtained earlier. Then we have 



+ a (los) 



^ = Y,B t (p 3 ,u q )d t v q + ^D(u q ,v q , Pj ) +( 2 (109) 



where & G ) and 6 G C^(]R b ) are 0(Z 9 ). 

By considering the differences between (106), (108) and (107), (109) one 
is lead to the following linear system of equations 

A°(U)d z (u - u q ) = A\U)d l (u - u q ) - (A°(U) - A°(u q ))d z u q 

+(A l (U) - A\u q ))d t u q + (G(U) - G(u q )) 



\c{U)-C{u q )+ jf s ( 



f{x?) r^ /rn , r ( o \ ,-i i 



(110) 
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6 

d z (v ~v q ) = Y J B\pj, U)di{v - v q ) + (B\ Pj , U) - B\ Pj ,u q ))d t v q 
i=l 

^(D(U,V, Pj ) - D(u q ,v q , Pj )) - 6 (111) 

for the unknowns v — v q ,u — u q with initial data v — v q = u — u q = 0. The 
quantities U and V are taken to be known smooth functions. Q q stands for 
the appropriate part of v q . 

4.4.5 Solving the field equations by a contraction mapping tech- 
nique 

By the theory of regular, linear symmetric hyperbolic PDE (Racke 1992) we 
can, given smooth (U, V) solve (111) for smooth v — v q and then (110) for 
smooth u — u q as long as U — u q , V — v q are 0(Z q ). These solutions exist as 
long as detg remains positive, where the elements g a/ 3 are extracted from the 
vector U. Equations (HO)-(lll) thus generate a map <E> : (U, V) — > (u,v). 
We will define a metric space (d, S) of functions in such a way that $ is 
a contraction mapping on S with respect to d. Some standard techniques 
from the theory of regular quasilinear symmetric hyperbolic systems (Racke 
1992) are used to achieve this. 

Let U defined on M 3 and V defined on M 6 belong to the space S iff the 
following hold: 

1. V(Z, x l , •) is supported in 0, for Z <T, where £1, relatively compact, is 
slightly larger than supp/°(x\ •) and supported outside B e (0) ioiZ < 
T with B e (0) slightly smaller than supp/ . 

2. The elements g af3 of U are such that \detg(Z, x { )\ > 5 VZ < T. 

3. u € c °°(o Q ), v e c^{o a x n). 

4. max{||?7 — u q \\ 8 , \\V — v q \\ s } < pZ q for Z < T, p const, s large. 

5. \\d z U\\s-i < L 

Here || || s is the L 2 type Sobolev norm of order s (Adams 1975). We 
choose s so large that the Sobolev imbedding theorem may be applied wher- 
ever needed in the sequel. 
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Lemma 4.2 If T is chosen sufficiently small and L sufficiently large 
then $ : S -> S. 

Proof First we deal with condition 4. Apply x % derivatives of order a to 
(110) and x l ,pj derivatives of order a' to (111) to get, after some rearrange- 
ment 

A°(U)d z V a (u - u q ) = A l (U)d l V a (u - u q ) + Ff - 

(112) 

d z V a '(v - v q ) = B i (x)d i V a '(v - v q ) + F?' - V a '6 (H3) 
where the Fj contain all the singular ^ terms as well as terms like 
-A^U^iA^U)- 1 A> (11)83(11 - u q )) 

etc. 

Now take the inner product of (112) with V Q (n — u q ) and of (113) 
with V a (v — v q ). Using the symmetry of A a ,B l and the uniform equiv- 
alence of the L 2 norm of V a (u — u q ) to 



/ 



A°(U)V a ( U - U q ) -V a ( U - Uq) (fx* 



in a standard way (Racke 1992) we sum derivatives of order < s to obtain 



||« " V q \\ s (Z) < C [ (\\V - V q \\ s + ^ \\ F ?h) (*) + ^ & 

J ° \a\<s (H4) 

\\u-U q \\ s (Z) < C [ (\\U-U q \\ s + \\F? || 2 ) (t) + t q dt 

Jo \*'\<s (115) 

where we also used ||^j|| s < cZ q . 

For the terms in F 2 a coming from B % one may readily use a Moser type 
estimate (Racke 1992) to obtain an I? bound in terms of the H s norm 
of v — v q , U — Uq. For the terms coming from D it is possible to obtain 
an L 2 bound in terms of the H s norm of U — u q , V — v q multiplied by ^. In 
this way (114) leads to 



\v ~ Vq\\ s < C 



(■z x 

/ Ik - V q \\ B + \\U ~ Uq\\ s + ~(\\U - Uq\\ s + \\V - V q \\ s ) + t q dt 

Jo f (116) 
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Similarly the terms in Ff coming from A a , G may be bounded in a standard 
way by ||u — Uq|| s , ||?7 — u q \\ s . Terms coming from C can, by Moser and 
Sobolev inequalities be bounded by a constant times Z~ l \\U — u q \\ s . The 
bound (iv) on U gives a bound for the speed of propagation in equation 
(111) and thus a bound on the support of Q(v) in the Pi variables. This 
entails that the momentum integral of Q may be estimated by \\v — v q \\ s , 
via the Cauchy-Schwarz inequality. In this way (115) leads to 

rZ 1 
\\u ~ Uq\\ s < C \\u- U q \\ s + \\U - U q \\ s + -(\\v - V q \\ s + \\U - U q \\ s ) + t q dt 

Jo t (117) 

Now write X = \\v — v q \\ s + \\u — u q \\ s and Y = \\U — u q \\ s + \\V — v q \\ s . 
Adding (116) and (117) gives 

X(Z)<C f X + Y + t q + -(X + Y) dt (118) 
Jo t 

From (111) we see that 



d n (v — v q 



= 

z=o 



for n < q. Thus v — v q is 0(Z q ). Similarly u — u q is 0{Z q ). 
Now put 

X\Z) = Z~ q X(Z) \\X\\' S = sup t- q \\X\\ s (t) (119) 

0<t<T 

with similar definitions for Y. Equation (118) now implies 

X\Z) < C [z' q X'(t q + ti' 1 + Y'{t q + t 9 - 1 ) dt + Z 



(120) 



and thus 



l*K^(iTi + ;)l*llKiTi + i 1 ' + T 



(121) 



By choosing q large and T small we can therefore arrange that \\X\\' S < p. 
It follows that $ preserves condition 4. Conditions 1, 2 and 3 follow from 
standard theory for linear hyperbolic PDE and from the finite speed of 
propagation inherent in the equations. 
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Now given condition 4 for (U, V) it is elementary, using a Moser estimate 
to show that 

\\d z {u - u,)||,_i < Ci (||« - uj s + (1 + Z~ l )\\v - v q \\ s ) + C 2 

(122) 

But by the preceeding result the right hand side can be bounded by a 
constant. Thus if L is chosen large enough then <& preserves condition 5. 
Hence $ : S -> S. 

Now define a distance function d on S according to 

d(p,p') = sup t- g \\p-p'\\ 2 (123) 

0<t<T 

Lemma 4.3 For small T and large q the mapping $ is a contraction 
on S with respect to d. 

Proof Let (U,V), (U',Y') G 5 and put = $(^7,F), («',*/) = 

$([/', V). Also write AU = U — U' etc. Then by a standard energy argu- 
ment one gets 



|Av|| 2 <C f \\Av\\ 2 + \\AU\\ 2 + 7(||AC/|| 2 + ||Ay|| 2 ) dt 
Jo t 



(124) 



and 



||Au|| 2 <c/ ||Au||2 + ||A^||2 + 7||Av|| 2 rft (125) 
Jo t 

Now put X = \\Av\\ 2 + ||Ait|| 2 , Y = \\AU\\ 2 + \\AV\\ 2 . Then 

X<c[ X + Y + -(X + Y)dt (126) 
Jo t 



If we now define 



\X\\' 2 = sup t- q \\u\\ 2 (127) 

0<t<T 



with a similar definition for Y then (126) leads to 

|Msc (?ri + i) M+c (?Tin) l|y|li (128) 

For small T and large q one may thus arrange that ||X|| 2 <(1 — <5)||y|| 2 for 
some 5 > 0. Hence the Lemma. 
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Now it is elementary to show that the set D S of functions p satisfying 
the condition 

SUp t _9 ||p||2 < p 

0<t<T 

is complete with respect to the metric d. Hence there exists p G S such 
that p = <&{p). Now let p k = ® k (p°) for some p° G S. Then d(p k ,p) — ► 
as /c — > oo and hence p fe is Cauchy with respect to d. 

By the Gagliardo-Nirenberg inequality (Racke 1992) one now obtains 

nz-v-p m )iU' < \\z- q {p k -p m )\\i ,s \\z~ q {p k -p m )\\\~ s ' /s 

(129) 

for s' < s. 

But $ : S — > S and so the ||Z -<? p fc || s are uniformly bounded. It follows that 
if we define the distance function d s > by 

d s 'ip,p') = sup t~ q \\p- p'\\ s i 

0<t<T 

then p k is Cauchy with respect to d s >. Hence \\p k — p\\' s , — > as k —>■ 
oo and \\p\\' s i < p. 

From the field equations one sees that the d z p k converge in C([0, T],H S _1 ) 
and therefore p G C([0, T], i/ s ') n ^([O, T], H 3 '' 1 ). By the Sobolev imbed- 
ding theorem p is classically C 1 and satifies the field equations ( p stands for 
the whole solution here). Since p is 0{Z q ) we can repeatedly differentiate 
the field equation to get p G f|p=o CP (.[°, T],H s '- p ). It follows that the parts 
of p corresponding to /, g a p and /i Q( g 7 belong to f]p =0 C p ([0,T], H s '~ p ). We 
could of course have chosen s' arbitrarily large and we thus get that /, g a p and 
ha/3-y belong to C([0, T s ], H s ), for T s < T, s > 0. The data at T s can then be 
used to construct data for the regular Einstein- Vlasov equations. Now it is 
well-known that the evolution of H s data for the regular EV equations stays 
in C([0, Z], H s ) as long as the C 1 norm remains bounded (Rendall 1997). 
In this way we see that our solution of the conformal field equations lies 
in C([0, T], H s ) for every s and is therefore C°°. 

It remains to show that the part Q a p of our solution agrees with the 
original definition (96)- (97). If we use the parts / and g a p of our solution 
to construct a quantity Q a p via (96)-(97) then a calculation gives that 

6 A 1 

dz(Qa(3 ~ Qaf3) = ^M l {u,Pj)di(Q al3 - Q af3 ) - —(Q a f3 ~ Qa?) 

i=i (130) 
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where the M l are symmetric. 

Since the singular part of this equation is 'negative definite' an L 2 energy 
estimate gives that Q a p = Q a p as required. 

Also note that for small T' we have F(x) = = 1 on the domain 

of dependence of 0' a x [0, T']. Thus u solves the conformal Einstein- Vlasov 
equations on 0' a x [0, T'\ and is determined by the right data. 

4.4.6 Uniqueness of solutions 

We now show that the solution p constructed above is the only smooth 
solution of the conformal field equations with the given data. 

Suppose then that we have two solutions {f,g a p), (/><7a/?) of the lo- 
calised equations (49)- (50) with the same Cq° data. These two solutions 
will then satisfy equations (96)-(97), (106)-(107). Let u stand for (g, dg) and 
let v stand for (/, Q) with similar definitions for u, -0. Since the system (106)- 
(107) is symmetric hyperbolic and the supports of / and / remain bounded 
on any closed time interval of existence we may obtain the following I? en- 
ergy estimates for the differences Au and Av: 

\\Au\\ 2 (Z) < C [ ||Au||2 + -(||Au|| 2 + ||At;||2)ds (131) 
Jo s 




Now letting X = ||Au|| 2 + ||Aw|| 2 we get that 

f z 

X(Z) <C (1 + s'^X ds (133) 
Jo 

We know from section 4.4.2 that the Taylor series of a solution to the field 
equations is determined by the initial data. Thus our two solutions (u, v) and 
(n, v) must agree up to arbitrary order in Z. This is to say that X < \ q Z q for 
arbitrary q. Plugging this into (133) gives 

X(Z) < C\ q Z« (1 + -^j (134) 

Now choose q > C and choose T small enough so that the right hand side 
of (134) is less than (1 — e)\ q Z q on [0, T] for some e > 0. We see iteratively 
that X < (1 - e) n X q Z q on [0,T]. Thus Au = Av = as required. 
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4.4.7 Dependence of solutions on initial data 



Suppose now that we have a 1-parameter family of small perturbations 
f®(x,p) of some initial datum /q depending smoothly on e,x and p for e £ 
(—1,1). We wish to show that the corresponding solutions (fe,gfj) of the 
conformal field equations depend smoothly on the parameter e. First it 
follows from an analogue of Theorem 3.1 that the family of initial met- 
rics h\j depends smoothly on e and x l . Now we may reconsider the field 
equations (136)-(141) for f,g,h,Q with e simply acting as an extra co- 
ordinate. By replacing the spaces H S (R 3 ) and H S (M®) respectively with 
H S (R 3 x (—1, 1)) and H S (E. 6 x (—1, 1)) we can now repeat the analysis of 
section 4.4.5 to get that the solution (/, g a p) depends smoothly on e. 

The results of section 4.4 may now be summarised as follows: 

Theorem 4.1 Given a smooth positive function f° on the cotangent 
bundle of a paracompact 3-manifold M satisfying the integral constraint 
(13) and such that for x in a compact set 

1. supp/°(x, •) C Sli, f° r some relatively compact Qi C M 3 

2. E 3 \supp/°(x, •) D $^2, for some neighbourhood of the origin C K 3 . 

there exists a unique solution (M x (0, T) , /, g a b) of the Einstein- Vlasov equa- 
tions with an isotropic singularity having f° as the orthogonal projection 
of the limiting particle distribution function at the singularity surface. The 
solution obtained depends smoothly on the choice of initial data. 

Appendix 

Theorem A.l. Consider a PDE of the form 



with A a symmetric and ^4° positive definite. 

If C{u)u° = and (^4°) _1 C(n°) has no positive integer eigenvalues then 
there exists at most one smooth solution u of (1) with u(0) = u°. 



A 



(u)d t u = A^ujdiU + B(u)u + -C(u)u 
u(0) = n° G C °°(]R n ) 



(1) 
(2) 
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Proof. Suppose we have two smooth solutions u, v of (1) with u(0) = 
v(0) = u°. One may now subtract the equation satisfied by u from that 
satisfied by v and obtain the following energy estimate: 

/"* „ 1 

\\u — v\\2(t) < K / \\u — v\\2-\ — \\u — v\\2 ds (3) 
Jo s 

on the time interval [0, T] for some constant K. 

Next note that the Taylor series of any smooth solution of (1), (2) is deter- 
mined by u°, by the following argument: 
Write u = u° + tu 1 + t 2 u 2 + ... + t q R(t, x l ). Then (1) implies 

+ 2tu 2 + . . . + qt q - l R + t q R'} 



-A\u)diU - B(u)u - {u 1 + tu 2 + t q ~ l R] = (4) 
Evaluating (4) at t = gives 

{A°(u°) - Ciu ))^ = A^u^diU + B{u°)u 
and thus u 1 is determined by u°. 

Now we differentiate (4) p times with respect to time and evaluate at t = to 
get an equation of the form 

(pA°(v°) - C(u°))u p = F{u p -\u p - 2 ...) 

and thus u p is determined inductively by n°. 

Since our two solutions u, v of (1) have the same Taylor series it follows that 
we must have the inequality ||n — v\\2 < \ q t q for t £ [0, T] and q arbitrary. 
Substituting this into the right hand side of (3) leads to 

\\u-v\\ 2 <K\ q t q (- + -^—) (5) 

Now choose q > K and choose T so small that the right hand side of (5) 
is less than (1 — e)\ q t q for some e > 0. Iteratively one gets \\u — v W2 < 
(1 - e) n \ q t q on [0,T]. Thus u = v as required. 
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